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Abstract 

We deal with the problem of reconstructing material coefficients from the far-fields they 
generate. By embedding small (single) inclusions to these media, located at points z in the 
support of these materials, and measuring the far-fields generated by these deformations 
we can extract the values of the total field (or the energies) generated by these media at 
the points z. The second step is to extract the values of the material coefficients from 
these internal values of the total field. The main difficulty in using internal fields is the 
treatment of their possible zeros. 

In this work, we propose to deform the medium using multiple (precisely double) and 
close inclusions instead of only single ones. By doing so, we derive from the asymptotic 
expansions of the far-fields the internal values of the Green function, in addition to the 
internal values of the total fields. This is possible because of the deformation of the 
medium with multiple inclusions which generates scattered fields due to the multiple 
scattering between these inclusions. Then, the values of the index of refraction can be 
extracted from the singularities of the Green function. Hence, we overcome the difficulties 
arising from the zeros of the internal fields. 

We test these arguments for the acoustic scattering by a refractive index in presence 
of inclusions modeled by the impedance type small obstacles. 
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1 Introduction and statement of the results 


1.1 Motivation of the problem 

Let n be a bounded and measurable function in such that the support of (n — 1) is a bounded 
domain We are concerned with the acoustic scattering problem 

(A + = 0 in (1.1) 


where := + Id* and Id* is an incident held satisfying (A + k^)V^ = 0 in For simplicity, 

we take incident plane waves V^{x,d) := e***"*^'® where 0 G and is the unit sphere. The 
scattered held Id* satishes the Sommerfeld radiation condition: 


dV^ 

d\x\ 


- mid* = o 



x\ —>■ oo. 


( 1 . 2 ) 


The scattering problem (1.1 1.2) is well posed, see 


Applying Green’s formula to Id* 


we can show that the scattered held ld*(a;, 6) has the following asymptotic expansion: 


„iK,\x\ 

V^(x,9) = - ——rld“(h, 6*) + 0(b|“^), b| —>■ oo, (1.3) 

4:71\x\ 

where the function ld“(a;, 6) for {x, 0) G x is the corresponding far-held pattern. 

Our interest in this work is related to the classical inverse scattering problem which consists 
of reconstructing n(x), x E fl, from the far-held data V°°{x, 9) for some {x, 6*) G x This 
type of problem is well studied and there are several algorithms to solve it in the case when 
X and 9 are taken in the whole see [26l [23, [28]. It is also known that this problem is very 
unstable. Precisely, the modulus of continuity is in general of logarithmic type, see PU] . 

Recently, based on a new type of experiments, a diherent approach was proposed, see for 
instance [SIEI and the references therein. It is divided into two steps. In the hrst one, we 
deform the acoustic medium by small inclusions located in the region containing the support of 
(n — 1), and measure the far-helds generated by these deformations. From these far-helds, we 
can extract the total helds due to the medium n in the interior of the support of (n — 1). In the 
second step, we reconstruct n{x) from these interior values of the total helds. Recently, there was 
an increase of interest in reconstructing media from internal measurements. This is also related 
to the hybrid methods introduced in the medical imaging community, see HI El El El EH ED 
for different setups and models. In contrast to the instability of the classical inverse scattering 
problem, the reconstruction from internal measurements is stable, see [HEiiiniEn]. However, 
the disadvantage with internal helds, i.e. internal values of the total helds, is the existence of 
their zeros which need to be properly dealt with to stabilize any algorithm. To overcome this 
short coming, one can think of using measurements related to multiple frequencies ki, 
see [2]- 

Our objective in this paper is to propose an alternative to overcome the last disadvantage. 
For this, we propose to deform the medium using multiple (precisely double) and close inclusions 
instead of just single inclusions. Then, measuring the generated far-helds by these multiple 
inclusions, we can extract, not only the internal total helds, but also the internal values of 
the Green’s function related to the non deformed medium n. This is possible because by 
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deforming the medium with multiple and close inclusions we generate scattered helds due to 
the multiple scattering between these inclusions. Precisely, the far-helds we measure encode at 
least the second order term in the Foldy-Lax approximation and not only the hrst order (or the 
Born approximation) as it is done when deforming with single or multiple but well separated 
inclusions. Finally, we extract the values n{x)^ x in from the singularities of this Green’s 
function. Hence, we avoid the problems coming from the zeros of the internal total helds. 

The accuracy of the reconstruction is related to the minimum distance d between the pairs 
of inclusions, in addition to the radius of the perturbations a. In practice, the perturbations are 
created using focusing waves, see ©[niEi] for instance, hence we are limited by the resolution 
of the used waves (acoustic waves for instance). In other words, the small perturbations cannot 
be very close. In our analysis, the ratio | is of the form with a parameter t G [0,1], and 
it can be chosen small by taking t near 0. This allows to have the minimum distance between 
the inclusions quite large compared to their sizes. Since the accuracy of the reconstruction, 
see Theorem 1.2, is of the order a*, with t in [0,t*] with some t* < 1, then one should hnd a 


reasonable balance between the limited resolution of the focusing waves, to be used to perturb 
the medium, and the desired accuracy of the reconstruction. More details related to this issue 


are provided in Remark 1.3 


We complete this section by describing the kind of deformations and the corresponding 
far-held measurements we use and then we state the derived formulas to extract the values of 
the index of refraction n. In Section 2, we justify these formulas with the emphasize on the 
explicit dependence of the error terms on the parameters modeling the small perturbations (i.e. 
the maximum radii, the minimum distance between them and the scaled surface impedances). 
In Section 3, we discuss the stability issue by providing the corresponding formulas when the 
measured data are contaminated with additive noises and when the small perturbations are 
shifted from their exact locations (as moving from step 2 to step 3 in collecting the data, see 
Section 1.3). We provide the regimes under which the derived formulas are still valid. 


1.2 Deformation by multiple inclusions 

We give the details of this approach by using inclusions of the form of obstacles of impedance 
type. Let B be the ball with the center at the origin and radius 1. We set Dm '■= clB + Zm 
to be the small bodies characterized by the parameter a > 0 and the locations Zm G 
m = 1,..., M. We denote by f/® the acoustic held scattered by the M small bodies Dm C 
due to the incident held t/* (mainly the plane incident waves U^{x,9) := V^{x,9) := 
with the incident direction 9 G where being the unit sphere), with impedance boundary 
conditions. Hence the total held := W + satishes the following exterior impedance 
problem for the acoustic waves 


(A -I- K n {x))U = 0 in U Dm i , 

V m=l 


M - 


dU* 

dlZm 

dW 

d\x\ 


+ 


= 0, 1 < m < M, 


dD„ 


inll^ = o 1— r ) fI cx). 


X 


(1.4) 

(1.5) 

( 1 . 6 ) 
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The scattering problem 
negative, see [T^ . 

Definition 1.1. We define 

1. d := min dmj, where dmj '■= dist{Dm, Dfi. 

^max OS the upper bound of the used wave numbers, i.e. n E [0, ^max]- 

The distribution of the scatterers is modeled as follows: 

3. the number M := M{a) := 0(a~^) < MmaxCi~^ with a given positive constant Mmax- 

4- the minimum distance d := d{a) ~ a*, i.e. dminof < d{a) < dmaxof, with given positive 
constants dmin cmd dmax- 

5. the surface impedance \m '■= , where Am,o 7^ 0 and might be a complex number. 

Here the real nnmbers s, t and fi are assnmed to be non negative. We call M^ax, dmin, dmax 
and Kmax the set of the a priori bounds. In ([12], Corollary 1.3), we have shown that there 
exist positive constants oq, A_, A+ [^depending only on the set of the a priori bounds and on 
rimax ■= ||n||Loo(Q) such that if 

a<ao, |A^,o| < A+, |3f?(A^,o)| > A_, fi<l, s<2-fi, ^<t, (1.7) 

then the far-held pattern U°°{x,6) has the following asymptotic expansion 

M 

6) = V^{x, e) + Y, y\zm, + 0{MCa), ( 1 . 8 ) 

m=l 

uniformly in x and 6 in Here C := rnaXmCm- The constant appearing in the estimate 0{.) 
depends only on the set of the a priori bounds, A_, A+ and on Umax- The coefficients Qm, 
m = 1, ...,M, are the solutions of the following linear algebraic system 

M 

Qm + CmG{Zm, Zj)Q.j = —CmV^{Zm,0), (1-9) 

j^m 


(1.4)-(1.6) is well posed, see [M [T3]. and we can also allow ^\m to be 


for m = 1,..., M, where 


C™, := 


Am I dDm I 

1 T Am^n 


( 1 . 10 ) 


and G(-, •) is the Green function corresponding to the scattering problem (1.1 1.2) 


^In the case where s < 2 — /3, we need no a priori conditions on A_, A+. 
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The quantity Im '■= Ji 


is equal to the radius of Dm, 
condition: 


dD^ 47r|s„-4| 


i.e. 


dsm, t € dDm, is a constant if Dm is a ball and this constant 
Im = a. The algebraic system (1.9) is invertible under the 


s < 2 - /3. 


( 1 . 11 ) 


We use surface impedance functions of the form Am = Ao.mO ^, i-e. with Ao,m 7 ^ 1 in the 
case /3 = 1 , so that the constants 


Cm, = 


dvr ^ A 


0,m 


— 1 + Aq ^ 


( 1 . 12 ) 


are well dehned. In the subsequent sections, we choose /3 = 1 and Ao,m := 1 — h > 0. In 
this case 

C'm =-47ra^“'^(l-a'*). (1.13) 


Hence (1.8) becomes 

17 “(£, e) = e) + Y. ^ ’ 


M 


(1.14) 


m=l 


uniformly in i. and 6 * in S^. 

1.3 The extraction formulas 

We proceed in three steps in collecting the measured data: 

1. First step. We measure the far-helds before making any deformation. In this case the 
collected data are V°°{—9i,9j), i,j = l,...,A^"o > 2. 

2 . Second step. We deform the medium by a isolated inclusion and then measure the 
corresponding data U°°{—9i, 9j) :=U°°{Zm; —9i, 9j), for the directions 9j, j = 1,..., Nq > 2 , 
keeping in mind that we redo the experiment by moving the inclusion centered at Zm for 
m = 1, 2,..., M and Z = 1, 2, see Fig|^ 

3. Third step. We deform the medium using two inclusions close to each other i.e. the two 
points Zm and Zm in fl, are such that \Zm — ^ml ~ as a —)■ 0, with f > 0, see Fig|^ In 
this case, the collected data are W°°(—9i,9j) ■=W°°{Zm; —9i,9j), i,j = 1 , ...,No > 2 . As 
in step 2, we create one couple of two close inclusions at once. Then, we move the couple 
of inclusions, for m = 1 ,..., M, inside fl. 


The reconstruction formulas of the index of refraction from the described measured data 
are summarized in the following theorem. 
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Figure 1: The figure describes the experiments in step 2. We create only one inclusion and 
measure the corresponding far-fields. Then, we move the inclusion inside hi. 



Figure 2: The figure describes the experiments in step 3. We create double inclusions distributed 
close to each other and we measure the corresponding far-fields. Then, we move the two 
inclusions inside hi. 
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Theorem 1.2. Assume n to be a measurable and bounded function in such that the support 
of {n — 1) is a bounded domain 12. Let the non negative parameters (3, h and t describing the 
set of the small inclusions, embedded in 12, be such that 

(3 = 1, 0 < h and h + 2t < 1. (1-15) 

Let m = 1,2,M be fixed. We have the following formulas: 

1. The vector of total fields (or can be recovered from 

as follows: 

C-] [u^{zL-e^,e,) - v^{-e„e,)] = -v\zl„e,)v\zl„e,) + o{a), (i.ie) 

fori, j = 1,2,..., No] I = 1,2. ^ As the two points z]^ and z^^ are close then the sign 
before the vectors above is the same for I = 1 and 1 = 2. 


2. Knowing the vectors of total fields for I = 1,2, (or for 

I = 1,2) the values of the Green function G{z(,,^, 2 ^) can be recovered as follows. We have 

Gizln,zK} := {G)i^2 , (1-17) 

where the matrix G is computable as 


VV^)-^C-^V\W^-Y^]V^{VV^)-^+C-^ = G+ O + a^) (1.18) 


with the matrix V constructed from (1.16) via {V)ij := V^{zl^,6j), I = 1,2 and j = 
1 , 2, ...,No and C := diagjCm ;}■ FI The matrices W“ and V“ are defined as (W“)j ,■ := 
W^{zl]-9,,e^) and (V“),,/:= fort,j = l,...,No. 


3. If, in addition, n is of class C", a G (0,1], therefrom the values of the Green function 
reconstruct the values of the index of refraction n{z\^) as follows: 



-Gizl 

IK 


lK\Zt 


+ 0{af), 2 > 0. 


(1.19) 


Remark 1.3. We make the following observations: 


1. We observe that due to the error 0{af), form = 1, .2, ...,M, we do not distinguish between 
n(z,^) andn^zffj since \zlj^ — z(^\ ~ as a —)■ 0 and hence \n{zff)—n{z‘ff)\ = 0{af) (in the 
best case a = 1). Hence these formulas provide the values of n{zm), for m = 1, .2,..., M, 
where Zm := z^^ for either Z = 1 or I = 2. 


^We recall that Cm.i = —47ra^“^(l — H), see (|l.l3). 

^ The matrix diag{Cm,z} is diagonal and the diagonal values are defined by the vector {Cm,i, Cm, 2 }- Since 
Cm.i = —47ra^“^(l —a^) as the shapes of the inclusions are the same, i.e. balls of radius a, we have diag{(7m,z} = 
—47ra^“^(l — a^)l 2 where I 2 is the identity matrix in 


^We can also take a = 0, i.e. n is only continuous, and obtain ( 1.19 ) with an appropriate change in the error 
term. 
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2. From the formulas (1.16), (1.18) and {1.19), we see that the error in reconstructing n 
is of the order hence we need h > 0 and 0 < t < The ratio 

between the radius a and the minimum distance d ^ is of the order | << 1 as 

a « 1. This rate becomes significantly smaller when t approaches 0. This means that the 
small perturbations do not need to be very close to each other to apply the reconstruction 

This makes sense and can be useful in practice since, due to 


formulas in Theorem 1^ 


the resolution limit, it is delicate to use focusing waves to create close small deformations. 
The price to pay is that we will have less accuracy in the reconstruction, see the formula 


(1.19) for instance. Hence, one needs to find a compromise between the desired accuracy 


of the reconstruction and the permissible closeness of the small perturbations created by 
focusing waves. For instance taking t = h = the error in the reconstruction is of the 

1 3 

order and the ratio ^ is of the order , as a « 1 . 

3. The importance of the parameter t, describing the closeness of the small inclusions, ap¬ 


pears in the last step (i.e. the point 3. in Theorem 1.2) to compute the values of n using 


the singularity analysis of G. Indeed, in (1.18), we can even take t = 0 with h G (0, 1). 
If we need to choose t very small (for practical or experimental reasons), then we can 
also use other methods to compute n{x) from G{x,z). For instance, from the Lippmann- 
Schwinger eguation G{x, z) = $(x, z) + ~ x)dt, x,z G 91, we can 

compute (nf — l)(x)G(x, z) by inverting an integral eguation of the first kind having as a 
kernel ^{x,z). Knowing G{x,z) and ipnf — l)[x)G{x, z) for a sample of points x,z &fl, 
we can reconstruct n from the formula {n^ — l)(x) = ■ The denominator 

G(-,z) is not vanishing if we choose appropriately different values of the source points z 
in 


4 . In this work, we use localized perturbations modeled by small inclusions with impedance 
boundary conditions. However, the natural localized perturbations should occur on the 
coefficients of the model the index of refraction n in our model. In this respect, 

using surface impedance type inclusions might be more attractive in the academic level 
than in practice unless we accept and use the following argument. In case where there is 
a large contrast between the medium and the background, and according to the Leontovich 
approximation, the contrast can be replaced by the surface impedance boundary condition. 
This means that we should use localized perturbations creating guite high contrasts. 


2 Proof of Theorem 11.2 


2.1 Step 1: Use of single inclusion to extract the field 

We create one single inclusion lo cate d at a selected point in This means that we take M = 
1 in the asymptotic expansion (1.8), i.e. U°°{z’^,x,9) — V°°{x,9) = V^{zl^,—x)Q + O (a^“^) 
and Q = —G„i,iV^{z\^i9). Hence 

t/“(4; -9., 9,) - V"“(-9..9,) = -C™,,l^‘(i.9i)l/‘(4,9,) + O . (2,1) 

Using the backscattered data TT^{z\^, —9j, 9j) and V°°{—9j, 9j), and since Gm,i ■= —47ra^“^(l — 














( 2 . 2 ) 


a") is known, from (2.1) we derive the formnla 


Hence we can constrnct either V^{zl^,9j) or —V^{zl^,9j) since the complex argnment will be 
compnted to an additive kn period. From (2.1), we see that we can compnte V^{zl^, 9i), i ^ j, 
with same additive period. Then, we can constrnct either the vector or the 

vector -{V\z'^,9j))f^^. 


Let ns now show why if zj^ is close to 
1 = 2 . 


2 


then this sign shonld be the same for I = 1 and 


We know that the total helds V*{-,9j),j := 1, are solntions to the Lippmann-Schwinger 
eqnations V^{x,9j) + f^{n‘^{z) — l)^{x, z)V{z,9j)dz = u^{x,9j), x G M^. Based on this 

eqnation and following the argnments in Lemma 2.2 of m, we show that V{x, 9) is in 
with the bonnd 




(2.3) 


where C depends only on the npper bonnd of n and f2, here CC 11. As 0 varies in the 
constant \\u\-,9)\\j^2(^q'^ is nniformly bonnded. We set C := (7snpggg 2 \\u\-,9)\\jj2(^Qy Hence we 
dednce that 

\v‘{zl,e)-v‘{zl,e)\ < c\zi-zl\. (2.4) 

Hence if the points z^^ and z"^ are chosen snch that ~ ^ml small enongh, then V^{zl^, 9) 
and V*{z^, 9) will have the same signj^ 


2.2 Step 2: Use of multiple inclusions to extract the Green’s func¬ 
tion e Q 

In this case, we create two inclnsions located close to each other, i.e. the two points z}^ and 
m = 1,..., M, in 11, are snch that \z‘^ — ^ml ~ a*, as a —)■ 0, with 1 > 0. Let m be hxed. 


From the asymptotic expansion (1.8), we have 


2 

W^{zl-X, 9) - V^{x, 9) = Y, V\zl, -x)QL + O {MCa ), 

1=1 

where 


(2.5) 








( 2 . 6 ) 


l' = l 
I'^l 


Hence, we rewrite (2.5) as 


W^{x, 9) - 0) = V ' (-h)B-W(0) + O {MCa ), 


(2.7) 


2.1 


^For m := 1, 2,..., M and F= 1, 2 fixed, the value V*{zl^, 9j) can vanish for some j but according to Lemma 
it cannot be true for ever j := 1,..., iVo where Nq is large. 
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where we use the vector V( 6 ') := 0), V*{z‘^, 0)]"'' and the 2x2 matrix B whose compo¬ 
nents are Bi^i/ := -G{zl^, z’^) for I ^ I' and := for / = 1 , 2 . 

We use a number N of directions of incidences (and hence of propagation directions) larger 
than the number of sampling points i.e. 2 directions, where we want to evaluate the 

index of refraction. With this number at hand, the matrix 


V:= [V(0i),V(02),...,V(0^)] 


( 2 . 8 ) 


has a full rank and hence VV is invertible. Precisely, we have the following lemma. 

Lemma 2.1. There exists iVo G N such that if the number N of incident directions 6j, j = 
is larger than Nq, i.e. N > Nq, then the vectors 'VhiOj) := [W(x^, 0^), W(x^, 
j = 1, ...,N are linearly independent. 


Proof. We observe, by the mixed reciprocity relations, that V^{zl^,9j) = G°°{—9j, zl^), where 




is the far-held in the direction —9j of the point source G{x,z\^ supported at the 


source point x. 


'^mJ 


is the Green function of the model 


We recall that the point source G{x, z: 

(1.1). With this relation at hand, the proof of the lemma follows the same arguments in 
where it is done for the vectors and j = 1, ...jN, with N large enough. We omit 

the details. Observe that V^{zl^,9j) = if the index of refraction n is equal to unity in 

□ 


Remark 2.2. In practice, one may take N = Nq = 2. This can be easily justified for the case 
when the contrast of the index of refraction, (n — l), is small. To see it we set, keeping the same 
notation as above, V := [V( 6 *i), V( 6 ' 2 )] whereN{9j) := [V^{zl,^,9),V^{z^,9j)]^, j = 1 , 2 . Corre¬ 


spondingly, we setY\{9j) : = 


ninz}„ -Oj ‘0-j 


T 


j = 1,2 andVi := [Vi( 6 *i), Vi( 6 * 2 ))]. A simple 


computation shows that the determinant o/Vi is e'‘ASi+S 2 )-zk,^^iK{z^-zh,)-d 2 _ ^iK.{z^-z^)-e 2 '^ ^ g 
if 6*1 7 ^ 6 * 2 . Now, we observe, by the Lippmann-Schwinger eguation for instance, that 
V = Vi-|-0(|n — 1|). Hence if the contrast n — 1 is small enough then the determinant 
ofN is not vanishing. 


Hence, by measuring the data and using the formulas 


(W~),,, := t,j = 1 ,...,Nq 

(2.9) 

and 

(V“),i:= {V^{-9,,9fi), t,j = l,...,NQ, 

( 2 . 10 ) 

we obtain the formula 

W°° - V°o = V’^B-^V + 0 (a^-^) . 

( 2 . 11 ) 

Observe that B^^ = — ^/ -|- CG j C where C := diag {Gm,i) and ^G j 

:= 0 , / = 1 , 2 and 

i 

(Gj := G{z\^,z''^,l fiz Ik Hence we can write 



B-i = -C + CGC + C*(max{G^,i}||CG^||2) 


( 2 . 12 ) 
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where ||||2 is the / 2 -norm for matrices. 

We recall that we have chosen the points and z^, in hi, such that \z^ — zl^\ ~ a* with 
t G (0,1), as a —)■ 0. With such a choice, we see that 


B-^ = -C + CGC + O 


(2.13) 


as 0(maXi{G^,z}||CG|||) = 0{a^ 

Finally, from the matrix B we derive the values of the Green function at the sample points 
^ = 1 , 2 , in G: 





(2.14) 


Indeed, combining (2.11) and (2.13), we deduce that 


W°o - V“ = -V^CV + V^CGCV + O (2.15) 


where V ' CV ~ and V ' CGCV ~ a 


,2-2h-t 


Remember that from step. 1, we compute 


±V (i.e. we do know the correct sign), but as in (2.15) we multiply the two sides by V and 
V^, then the missed sign of V is irrelevant. 


The approximation (2.15) will make sense, i.e. the error term is dominated by —V^CV + 
V^CGCV, a 3 — 3h — 2t > 2 — 2h — t and 2 — h>2 — 2h — t, i.e. 0 < t + h < 1. In add ition , 
V''^CGCV = 0(V''^CV) a 2 — 2h — t > 1 — h, i.e. t + h<l. Hence the approximation (2.15) 
makes sense if0<f + /;,<l. This last condition is satisfied since we assumed that 0 < h and 
h + 2t < 1. 


From Step 1, we showed already how to extract the matrices ±V, hence we can compute 
W“ — V°° — V''^CV . Since the matrix C is diagonal, with non zero entries, and the matrix 
VV^ is invertible, as V is full rank matrix, then we can recover the matrix G as 


C-i(VV^)-^V [W°° - V“] V^(VV^)-iC-^ + C-^ = G + 0(0^'^-^) O 

(2.16) 

or 

C-i(VV^)-^V[W“-V“]V^(VV^)-^C-^ + C-^ = G+ O + a^) . (2.17) 


We see that = o(G) ii h > —f,i.e. f + h. > 0 and a 


l-h-2t 


= o(G) li 1 — h — 2t > —t, i.e. 


t + h < 1. Hence (2.17) makes sense if 0 < f + /i < 1. As mentioned above, this last condition 
is satished since we assumed that 0 < h and h + 2t < 1. 

Finally, we set 

•= (G)i,2, (2-18) 


as the reconstructed values of the Green function G at the sampling points (z,(^,z^), m = 

1,2,...,M. 


2.3 Step 3: Extraction of the index of refraction from the Green 
function 

We start with the following lemma 
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(2.19) 


Lemma 2.3. We have the asymptotic expansion: 


G{x, Zj) 


^iKn{zj)\x—Zj\ 

47r|a; — Zj\ 


+ 0{\x 


Zj\), as \x — Zj\ —)■ 0. 


Proof. We know that (A + K^n^{x))G = —5, in and ^j{x,z) := ^ ' satisfies (A + 

K^n^{zj))^ = —S, in where both G and satisfy the Sommerfeld radiation condition. 
Then Hj{x,z) := {G — ^j){x,z) satisfies the same radiation condition and 

(A + K^n‘^{zj))Hj = K^{n^{zj) — n‘^{x))G, in (2.20) 


Mnltiplying both sides of (2.20) by <hj and integrating over a bonnded and smooth domain B 
we obtain: 


Hj{x,z) = / {n‘^{zj) — n‘^{t))^j{t,x)G{t, z)dt 


( 2 . 21 ) 


'B 


I dB 


Hj{t,z)d^(t)^j{t,x)ds{t) + / d^(t)Hj{t,z)^j{t,x)ds{t), 


'dB 


and then 


'VxHj{x,Zj) = —K^ / {n'^{zj) — n‘^{t))^j{t,x)'VxG(t, Zj)dt 


( 2 . 22 ) 


IB 


'dB 


Hj{t, Zj)Vxdu(^t)^j{t, x)ds{t) + / d^^t)Hj{t, Zj)Vx^j{t, x)ds{t). 


'dB 


Here, we choose B snch that CC B. Since {n?{zj) —n‘^(t))'VG{t,Zj) = 0{\t — Zj 


-2+a' 


,t e B, 


as the singnlarity of Va;G(f, —Zj) is of the order |f —and n^{zj)—n^{t) = 0(|f —Zj|“), since 
n is of class C", we dednce that J^{n^{zj) — 77,^(f))<hj(f, x)VxG{t, Zj)dt = 0{J^ \t — — 

x\~^dt) = 0(1), X & B since a > 0, see ([31], Lemma 4.1), for instance, for the snch integrals 
o f of singnlar fnnctions. The two last integrals appearing in (2.22) are of the order 0(1) for 
X eQ CC B. 

Now, let i? 2 j,r be the ball of center zj and radins r, r << 1. We use (2.21) applied to B^.^r 
instead of B: 


Hj{x,z) = / {n‘^{zj) — n'^{t))^j{t,x)G{t, z)dt 


(2.23) 




- Hj{t,z)d„(^t)^j{t,x)ds{t)+ dy(t)Hj{t,z)^j{t,x)ds{t) 

ddBz.,r ddBzj,r 

for X G B^.^r- We see that 

f {n‘^{zj) — n‘^{t))^{t,x)G{t, Zj)dt = C)(^ f \t — Zj\~^~^°‘\t — x\~^ds{t)^ , ior x G B^j^r 

*' ^Zj^v ^Zj,r 

which we can rewrite as 

{n‘^{zj) — n‘^{t))^{t,x)G{t, Zj)dt = o(^{ \t — Zj\~‘^~^‘^°‘ds{t))^ { j' \t — x\~‘^ds{t))^'^ , 
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for X G Hence 


{n^{zj) — x)G{t, Zj)dt = 0(r^+"), for a; G 


'B, 


Since VHj{t, Zj) is bounded for t E Q, then 

d^(t)Hj{t,Zj)^j{t,x)ds{t) = o(^ 


\^j{x,t)\ds{t)] = 0(r) 






for X E Hence, we can write (2.23) as 

Hj{x,Zj)+ [ Hj(t,Zj)d^(^t)^j(t,x)ds(t) = 0{r), for X E (2.24) 


IdB^ 


Taking the trace in (2.24) to dB^.^j. and using the trace of the double layer potential 
^ Hj{t, Zj)d^(^t)^j{t, x)ds(t) = Kr{Hj{-,Zj)){x) - \Hj{x,Zj), x E dB^.^r, where Kr de¬ 
notes the double layer operator, we derive the equation 


Zj)) = 0(r), for X E dB^.^r- 


(2.25) 


Then 


\\Hj{-,Zj)\\L 2 (dB,^,r) - O ^ ||(- + Kr) ^\\ciL2iaB,^,r), L2idB,.,r))^ ■ 

We have the following estimate from HBl: 

11(9 + tlli:(L2(aB9,,.), L2(as.,,.)) < IK^ + ^'■’i)^tlli:(i.20s}. L2(aB)) 


where B := Bo,i- Hence 


\Hj{-,Zj)\\L 2 idB..r) = 0{r). 


(2.26) 


Going back to (2.24), we estimate 
\Hj{x,Zj)\< j 

and then 


dB^,, 


\Hj{t,Zj)\\d^,^t)^j{x,t)\ds{t) + 0{r) 


\Hj{x,Zj)\ < ||-Hj(-,Zj)||L 2 (gB^ )|||5j9(p<hj(-,t)||i2(aB^ ) + 0{r). 


From (2.26) and the fact that ||5^(t)*hj(-,t)|| 1 , 2 ( 95 ,.,,) = \t - Zj\ ‘^dt)^ = 0(1), we 

derive the estimate 

Hj{x,Zj) = 0{r) 

for X G Bz r and in particular for x such that |x — Zj\ = r. Hence, from the dehnition of 
H, (x,j), we deduce that 

^iKn{zj)\x—Zj I 

G(x, Zj) = ——I-- -|- 0(|x — Zj\), as \x — Zj \ —>■ 0. 


47r|x — Zj 


□ 
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From (2.19), we deduce that, 


G(x, 


47r|a; — 


+ i K n{zm) +0{\x — Zm\), as \x — Zm\ ^ 0 (2.27) 


and then 


n(zrr,) = lim 


• G(x, ZjYi) . 

ZK ZK A7r\x — z„ 


(2.28) 


are the sampling points added to the small inclusions. 


The formula (2.28) can be used for x := z"^ and Zm ■= where zl^ and z^, m = 1,2,..., M 

+ 0(a*), f > 0. 


'^(Zm) 


G(Zjr^.i ^m) 


ZK 


m47r|z^-z^,j 


(2.29) 


3 The stability issue 


The object of this section is to analyze how the formulas are accurate when the measured far- 
helds are contaminated with noise and when the locations of the small perturbations z\^,m := 
1, ...,M and / = 1,2, are slightly shifted when passing from step 2 to step 3 in collecting the 
measured far-helds. Precisely, we will derive the regimes on the amplitude of the noise and the 
length of the shifts under which our formulas are still applicable. 


We recall the main formulas (2.17), (2.18) and (2.29); 


C-i(VV^)-^V [W“ - V'^] V^(VV^)-^C-^ + C-^ = G + O + a!^) 


and 


'^(Zm) 

Combining them, we obtain: 


G{Zra, Z^) (G)i,2, 

_ ^ _ 

ZK ^ ^^47^\zl-z}^ 


+ 0(a‘), t > 0. 


(3.1) 

(3.2) 

(3.3) 


iKn(zm) = [C-^(VV^)~^V[W“ - V“] V^(VV^)-^C-^+ - 

+0 + a^)+0 (a*) . 


- Z. 


n (3-4) 


We recall also that V = [V(0i), V(02), • • • ,V(0^)] where V(0) := [V\zi,9),V\zi,9)]^ 
and V*(zl^,9j) can be computed, see ( 2.1[ ) and (2.2), as 

C-][U^i^L■,-^^,9J)-V^{-9,,9,)] = -V\zl„9,)V\zl,9,)+ 0(a). (3.5) 


We base our remaining analysis on the two formulas (3.4) and (3.5). 

We set z)^ to be the possibly shifted positions of the small perturbations while moving from 
the second step to third step of collecting the far-helds, see Section |1.3[ We set also 
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r] ■= max \zi, - zl,] 


and 


d := 


mm 




where -D^’s are the ball of center 5^’s and radins a respectively. Let now 5^, 5u and 5^ be 
the additive noises corresponding to V°°, U°° and W°° respectively. The actnal measnred data 
is 14 °°, and satisfying respectively 


|rr(-9i,«j) - < iv, \ur{4n,-0i,ej) - < s. 


and 


\wrisL; -e.,e,) - -e.,e,)\ < h 


00 / 


(3,6) 


nniformly in terms of 6 i, dj 


i,j = l,...,A^"o. Observe that since 


C-; p“(4.i-9..9i) - = -V(4,9.)V(4„9 j)+ 0(a) then 

C-; [c/r(i. -»i, « 3 ) - 44(-»h Oi)] = -r'(4.. 9.)V(4. 4 ) + o (n) + + «„), ( 3 , 7 ) 

We set V = [V( 9 i),V( 92 ),," ,V(9j,)] where V(9) := [l/‘(4„,9),l/‘(5^,9)]^. Remember 
now that in the second step, the data are collected at the shifted points 5* , m — 1,2,A/ and 


/ = 1,2. Hence instead of (3.4), we have 


lK,n[Zm) = 


C-\YY )-^V 


W" 






V’(VV’)-iC'^ + C 


3-1 


4ir|4 


dl(3,8) 


+0 ( a, 


,l-h-2i 


+ aM + O a 


where we set W“ to be the matrix dehned as := W°°{zl^] —6i, Oj) and correspondingly 

we set the matrix W“ dehned as (W“)jj := —^j, 6j) . Here i satishes the conditions 


0 < h, 0 < f, h + 2f < 1. 


(3.9) 


From this formula, we get the corresponding one with noisy data W“ and 


IKniZm) = 


c-i(vv )-^v 


Wf -Yf 


ZT 


V (VV )“'C-^ + C 


3-1 


1,2 


^Az'L - zl 


7(3.10) 


+((5,, + (5,)|C-^(VV^)-iV||V(VV^)-^C-^| + O + +0 fa‘ 


Now, we need to replace the values of V by the values of V, more precisely the computed 
values of V from step 1 with noisy data. Let us set V 5 the computed one from step 1. We 


deduce then from (2.4) and (3.7) that 


|V 5 - V| = 0 {a + a'^-\6u + 5.) + n). 


(3.11) 
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With this estimate at hand, the hnal formula is: 


IKniZm.) = 




Wf - 


Vj(V,V^ 


+ c 


-1 


1,2 47r|5^ - 5^1 

(3.12) 


+0{a^^-‘^)0{a + + 5.) + il)\^T - V^l + (5,. + 5,)\C-\\\ ' )-'V||V(VV^)-'C 


-1/-I-1I 


+0 (a 


Z-h-2t 


+ ) + O [a 


As |W“ 


. V“| < |W“ 


■ W 


the points we have |W'^ 


+ 5v + 5w and from (2.15), applied to the points z^^’s instead of 
= 0(a^“^), then 




|Wr-Vri = 0(a'-") + 5^ + 5.. 

T, 


(3.13) 


\—I/"' —1 I 


or 


What is left is the estimate of the term |C“ (VV )“^V| (respectively |V(VV )“^C 
the term |(VV )~^ V| (r espectively |V(VV )“^| ) since the matrix C = —47ra^“^(l — a^)l 2 - 
Arguing as in Lemma 2.1, the matrix VV is invertible. However, what we need here is its lower 


bound. Wet set its lower bound. Then |C ^(Vv"'") ^V| = 0{a^ ^ = |V(VV') 

With these estimates, the error term in (3.12) becomes 




— I/"' —11 


0{a^^-^)0{a+a'^-\6u+6,)+r]){0{a^-^)+6^+6,)+{6^+6,)0{a^'^-y-^+ O a 


l-h-2t 


+ aM+0 a 


or 


0(a^)+0(a^^ ^)(5u+5^)+0(a^^ ^)((5u,+(5^)+0(a^^ ^)((5M+(5,(,)((5i„+5^)+0(a^ ^)r7+0(a^^ ‘^)Ti{6m-\-5v) 

{6^ + (5,)0(a2^-2)/-2 + O + O (a*) . 

If we reasonably assume that we have uniform noise for our different measurements, i.e. 

= 0{6),6y = 0(5) and 6^, = 0(5), then this reduces to 

0(0^*) + + 0{a^^-^)6^ + 0{a^-^)r] + 0{a^’^-^)r]6 + O + O (a') . 

We summarize these hndings in the following theorem 

Theorem 3.1. Assume n to be a measurable and bounded function in such that the support 
of {n — 1) is a bounded domain H. Let the non negative parameters (5, h and t describing the 
set of the small inclusions, embedded in be such that 

/3 = 1, 0 < h and h + 2t < 1. (3.14) 

We assume that the measured far-fields are contaminated with an additive noise of amplitude 
S, i.e. the actual measured data is and satisfying respectively 

\Vs^{-e^,e,)-v-^{-e,,e,)\<5, \ur{zl--e,,e,)-u-^{zl,-e,,e,)\<5 

and 

<« (3.15) 

uniformly in terms of 9i,9j G i,j = 1,..., Vq. 
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1. Let m = 1, 2,M he fixed. The vectors ±(^^*( 2 ;^, 0^), 9j)) satisfy the formulas 


C-', K(4. -9i, Oi) - VS(-9., 9;)] 


-V'‘(4,«i)V'‘(4„ 4) + O (o) + 0(o'>-‘)i. (3.16) 


2. Let -LVs be the constructed fields from step 1. We set rj to be the maximum shift of the 
centers of the perturbations from step 1 to step 2, that we denote by zl^,m = 1, ...,M and 
/ = 1,2. Let also the minimum distance between these shifted perturbations, i.e. d, he 
of the order a*, t > 0. Using the noisy data and we have the reconstruction 
formulas, I = 1,2, 




C-\V5Vl)-^Vs 


+0{o!^) + + 0 {a 


Wf -Yf 

3h-3)^2 ^ 


vJiVsV])-^cr^ + cr^ 

^-^)r] + 0{a^^-^)r]6 + O ( 


1 

47r|z^-z^| 

(3.17) 

^i-h- 2 i^ + o . 


The formula (3.17) will make sense if the error term is of the order o(l) as a —)■ 0. It is clear 
that the noise level shonld decrease in terms of a. We set 6 = 0(a'^^) and p = 0(a‘^^), then this 
error has the form 


O(a^) + 0(a'?i+2'*-2)C2 + ^ + O (aM . 


Hence the reconstrnction formnlas apply if we have the following relations: 


(3.18) 


0<h<l, 0<f< 


1 -h 
2 


qi > 2 — 2h, q2 > I — h 


(3.19) 


even if the constant ^ 
apply if 


might be qnite large. The formnlas (3.16), to compnte the total helds. 


0 < h < 1, gi > 1 — h. 


(3.20) 


We observe that the conditions on the noise level to reconstrnct the index of refraction (i.e. nsing 
both single and mnltiple pertnrbations) are more restrictive than those needed to constrnct only 
the total helds (i.e. nsing only single pertnrbations as it is snggested in the previons literatnre), 
which is natnral of conrse. 


4 Conclusion and comments 

We have shown that by deforming a medinm with mnltiple and close inclnsions, we can extract 
the internal valnes of the coefficients modeling the medinm from the far-held data. We have 
tested this idea for the aconstic medinm modeled by an index of refraction where the small 
inclusions are impedance type small scatterers. 

1. Let us consider the more general model V ■ c^V -|- where c is the velocity which we 
assume to be smooth and eventually diherent from a constant and the small perturbation 
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are also of impedance type R_ In this case, arguing as in [12] , we can derive the same 

asymptotic expansion as in (1.8)-(1.9)-(1.10). Hence proceeding as in the three steps we 
described above, we obtain: 

C-i(VV^)-^V[W~-V“]V^(VV’^)-^C-^ + C~^ = G+ O + a^) (4.1) 

and instead of (2.19|), we obtain 


G{x^ Zra) 


A'KC^{Zm)\x - Z„ 


+ Oi\n{\x - Zm\)), as |x - ZrnI-t 0. (4.2) 


From (4.2), we derive 


G(x, ZfYi'j 


_47rc‘^{zm)\x - z„ 


+ 0{\n{\x - Zm\)), as\x-Zm\^0- (4.3) 


Combining (4.3) and (4.1), with x = z^ and Zm ■= z^, we derive the formula 
c-^{zJ = A7r\zl-zl\ [C-\VV^)-W[W^ + (4.4) 

+0(a* In a) + + a’^) 

for z = 1,..., M, recalling that \zl^ ~ ^m\ — ^ — 1) M. 

We see that we can derive a direct formula linking the measured data W°° and V°° to 
the value of the velocity c on sampling points Zm, m = 1, 2,..., M in hi. The corresponding 
formulas using the noisy data can also be derived. 

2. The more realistic model would be V ■ c^V + where the small perturbations occur 
in the coefficients c and n and not as small impenetrable obstacles as we dealt with in 
the present paper. One needs hrst to derive the corresponding asymptotic expansion 
taking into account the parameters modeling these small inclusions and then use them 
to design a reconstruction method to extract c or/and n from similar measured data as 
we did in the present work. In the model we dealt with in this paper, we have chosen 
the perturbations to be spherical and the surface impedances to be of the form Am := 
(1 — a^)a~^. With this choice, the error in the approximation (1.14) allows us to recover 


the fields created at the second order scattering (and not only the Born approximation). 
This is the key point in order to reconstruct the index of refraction n using multiple 
perturbations. To extend this argument to the case when the small perturbations occur 
on the coefficients c and n, one needs to derive the Foldy-Lax approximation for such a 
model allowing to recover helds created at least at the second order (and, again, not only 
the Born approximation). This issue is related to the full justihcation of the Foldy-Lax 
approximation (i.e. at higher orders). The justihcation of the accuracy and the optimal 
error estimate in the approximation of the scattered held by the Foldy-Lax held is in its 


®As already mentioned in the point 4 of Remark 1.3 this is possible only if the localized perturbations create 
high contrasts so that, by the Leontovich approximation, we replace the contrast by impedance type boundary 
conditions. 
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generality a largely open issne bnt there is an increase of interest to nnderstand it, see 
for instance dni [m [131 El ESI [21]. 

Finally, we do believe that this approach can be applied to other mnltiwave imaging 
modalities as well. 
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